Kondo Insulator: p-wave Bose Condensate of Excitons 
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In the Anderson lattice model for a mixed-valent system, the d — f hybridization can possess a p- 
wave symmetry. The strongly-correlated insulating phase in the mean-field approximation is shown 
to be a p-wave Bose condensate of excitons with a spontaneous lattice deformation. We study the 
equilibrium and linear response properties across the insulator-metal transition. Our theory sup- 
ports the empirical correlation between the lattice deformation and the magnetic susceptibility and 
predicts measurable ultrasonic and high-frequency phonon behavior in mixed-valent semiconductors. 

PACS: 71.28.+d, 71.27. +a, 43.35.+d, 71.10.Fd 
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Mixed valent compounds containing /-electrons in 
their insulating phase possess very distinctive electronic 
properties showing strong electron correlation and 
are commonly called Kondo insulators, even though they 
are often not in the Kondo regime 0. The insulat- 
ing state may be regarded as a condensate of the slave 
bosons which are used to simulate the high energy cost of 
too many /-electrons on the same ion site [0. Through 
the /-level and conduction band hybridization, this Bose 
condensate drives the condensation of /-hole-conduction- 
electron excitons, creating a band gap. This theory has 
been used to explain a number of thermal, transport, and 
optical properties 0. 

In this paper, we adopt the same slave-boson treat- 
ment of the Anderson Lattice Model to explore the con- 
sequences of the d-character of the conduction electrons 
so that the d-f hybridization at wave-vector k exhibits 
inversion antisymmetry, i.e. 



V(-k) = -V(k) = V*(k). 



(1) 



This leads to odd-symmetry pairing of d-electrons and 
/-holes. (See below on how a p-wave symmetry can 
arise.) Coupling of the d-f excitons to the lattice defor- 
mation leads to a structural phase transition known as 
the ferroelastic transition ||. We deduce the tempera- 
ture dependence of the equilibrium properties (lattice dis- 
tortion and phase diagrams) and of the linear responses 
(including the magnetic susceptibility, elastic constants, 
phonon frequency and damping). These properties show 
the characteristic dependence on the p-wave condensate 
and their observation can be used to distinguish the s- 
or p-wave nature of the Kondo insulator state. We shall 
make the case below that some of the p-wave characteris- 
tics have already been found in certain Kondo insulators. 
The establishment of the p-wave characteristics in a class 
of Kondo insulators could be of importance in providing a 
central framework for understanding their unusual prop- 
erties in terms of the symmetry of the condensate. 

The contrast between the p-wave condensate in the An- 
derson lattice and the s-wave condensate in the Falicov- 



Kimball model J7j should also be noted. In the latter, 
the d-f hybridization is neglected in comparison with the 
strong d-f Coulomb attraction which is responsible for 
the exciton condensation. Since each d-f exciton carries 
a dipole moment, the condensate behaves like an "elec- 
tronic" ferroelectric. On the other hand, from symmetry 
considerations, the p-wave condensate does not produce 
a macroscopic electric polarization but can yield a fi- 
nite lattice distortion. Responses to lattice vibrations in- 
duce coupling to the collective modes of the p-condensate, 
which are fundamentally different from the electronic fer- 
roelectrics. When both the hybridization and d-f inter- 
action are present, the properties of the s and p wave 
coexist. Details will be presented in a long paper while 
here we concentrate on the limiting case when hybridiza- 
tion dominates, in contrast to the other limit H. 

We consider a model with one /-electron and one con- 
duction electron per unit cell p|. For simplicity, we 
include only the one component each of the d and /- 
bands which yield the hybridization obeying Eq. (Q) and 
neglect the / orbital degeneracy. To consider the for- 
mation of the condensate and the related dynamic re- 
sponse, we replace the onsite interaction in the Anderson 
model by the slave-boson contraints and further retain 
only zero-momentum slave bosons which are present in 
macroscopic numbers [0. The Anderson Hamiltonian is 
thus reduced to: 

H = J2 Kk) - n] 4A CT + (e f - A - n) /L/k. 

k(7 kcr 

+ -±= £Wk) dtj^tf + h.c.} + A {N s - fet 6 ) , (2 ) 
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where d\^ a , /k CT , b are annihilators for the conduction and 
/ electrons of wave vector k and spin a = ±1, and slave- 
boson of zero momentum. The conduction band energy 
is denoted by e(k) and the chemical potential by p. A is 
the Lagrange multiplier which enforces the constraint of 
at most one /-electron per cell for N s cells (on average). 

The equations of motion for the d-f exciton condensate 
are physically transparent upon expressing the density 
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matrix in terms of a pseudo-spin vector S^ 
S£ + iS£ = 2{dU klT ) 



(3) 



In the mean-field approximation, the spin vector satisfies 
a Bloch equation [Q: Sk = J^k x Sk, where the precession 
frequency is given by 
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F*(k) 



(b) 



F*(k) 



f(k) , (4) 



where 3? and 9 denote the real and imaginary parts and 
£(k) = \ [e(k) — e/ + A] . The expectation value at a finite 
temperature of the slave-boson condensate (b) is driven 
by the exciton condensate || (d kcr f klJ )- 



i5 t {6) + A(6) = -^^y(k)(4 <r / ktr ). 



(5) 



The zeroth-order (in the absence of external perturba- 
tion) solution of the above equations yield the thermal 
equilibrium state. The spin vector S^ ' is a unit vector 
in the direction opposite to the zeroth order precession 
frequency fl^ , yielding the exciton order parameter: 

(4JuJ 0) = aV*(k) lf(E k+ ) f(E k _)} /2£ k , (6) 



where f(E) is the Fermi distribution, a 
and the quasi-particle energies are Ek± = 



WfZ A<°) + 

e(k)] ± E k , where £ k = y/[^ 0) (k)] 2 + a 2 |F(k)| 2 . From 
Eq. (||), the exciton order parameter being driven by the 
hybridization V(k) must have the same symmetry. By 
Eq. ([!]), the pairing state is antisymmetric under inver- 
sion. More specifically, if V(k) is one of the p components 
of the hybridization matrix elements between the d and 
f states, the exciton order parameter would indeed have 
a p-wave symmetry. 

To evaluate the mean field phase diagram from the ze- 
roth order equations, we further adopt a simple model 
with the conduction band having a constant density of 
states with a bandwidth of 2W. The hybridization is 
approximated by V(k) = i Vq sgn(fc z ), where Vb is a real 
constant. The vanishing gap case is avoided by choos- 
ing the energy surface £(°)(k) = not to intersect the 
plane where V(k) vanishes. Some of the results are: (1) 
At T = for arbitrary €f, the ground state is always a 
p-wave condensate of excitons (correlated insulator). (2) 
For tf > 0, there is no finite-temperature phase tran- 
sition. (3) For ef < 0, there is a second-order phase 
transition at finite temperature between the correlated 
insulator and a simple metal. 
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FIG. 1. Condensate strength and spontaneous lat- 
tice deformation versus temperature for Ho = — 0A7W, 
N s = 3 x 10 27 /m 3 , C = 10 10 J/m 3 for the simple model whose 
other parameters are specified in the text. Squares are exper- 
imental results from Ref. 10 for Ce3Bi4Pt3. 

Figure 1 shows the temperature dependence of the 
slave boson condensate a(T) for e/ = —0.2 IF and Vo = 
0.3 W . The resultant fcgT c = 0.123 IF indicates this case 
to be in the intermediate coupling region. For the rest 
of this paper, the numerical results are reported for the 
same set of parameters. 

Two immediate consequences of the odd symmetry ex- 
citon condensate are the absence of polarization and the 
possibility of lattice distortion. For a more quantitative 
study, we include in the Anderson model a coupling with 
the strain field 9 through the deformation potential S(k) 
H by replacing V(k) by V(k) + H(k) 9. From the form 
of the electron-phonon interaction we have included, the 
deformation potential S(k) arises from the change of the 
hybridization V(k) to first order in the lattice displace- 
ment and, therefore, must have the same odd symmetry 
under inversion. The lattice deformation creates a stress: 



s = C9+- 7 =^Z(WU^)(b f ), 



(7) 



ka 



where C is the bare elastic constant. The zeroth order 
terms show that the exciton condensate induces a spon- 
taneous lattice deformation 9^°\ The feedback of the 
lattice distortion should be included in Eqs. ([|) and (||). 
A tight- binding derivation shows that 3(k) is the same 
order of magnitude as V(k) and the distortion 9^ is al- 
ways less than a percent (see below). Therefore, we may 
neglect the self-consistent inclusion of the lattice distor- 
tion in the determination of the exciton and slave-boson 
condensates. 

Figure 1 also shows the temperature dependence of the 
lattice distortion 9^°' in comparison with experiment on 
Ce3Bi4Pt3, chosen as an example of a class of small-gap 
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correlated semiconductors 10 IJ. The deformation po- 
tential is approximated in a similar way to the hybridiza- 
tion: S(k) = i So sgn(/c z ) where 2o is a real constant. To 
fit the experimental curve, we have chosen the param- 
eters given for our simple model and used T c = 350 K, 
which yields the conduction bandwidth of 2W ~ 0.49 eV. 
The overall agreement of the calculated curve with the 
measured curve should not be taken too literally but the 
support of the spontaneous structural phase transition 
accompanying the exciton condensation in the mixed- 
valent system in our theory is noted. This structural 
transformation belongs to a well-known class of ferroe- 
lastic phase transitions || . 

By adding a perturbing term to the equations of mo- 
tion, we can calculate the first-order changes. To calcu- 
late the Pauli spin susceptibility, we use the perturbation 



H B = -LIB J2^ d<7d l^a + 9fVfLf^a)B, 



(8) 
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where B is the magnetic field in a fixed direction, fis is 
the Bohr magneton, and gf and gd are phenomenological 
gyromagnetic ratios for / and d electrons. The Pauli 
susceptibility is 

X = Mb 2^ \ 2aV ( k )\ \9f-9d) (2E k ) 3 

+ (2/4/fc fl T) £ g 2 J(E ks )[l-f(E ks )}, (9) 

k,s=± 

where g± = gdA T (k) + gfA±(k) is the effective gyro- 
magnetic ratio for quasiparticles in the mixed (±) bands 
with spectral weights A±(k) = \{1 =p £ (0) (k)/£ k }- The 
first sum at the right-hand-side of Eq. (||) is due to inter- 
band mixing proportional to a 2 of the slave-boson con- 
densate, while the second sum is from the individual band 
contributions. Figure 2 and its inset show the temper- 
ature dependence of the spin susceptibility ^ and its 
relation with the lattice deformation as well as compar- 
isons with the experimental results for Ce3Bi4Pt3 |Io|,|l|. 
The observed low-temperature upturn of the suscepti- 
bility was sample dependent and was ascribed to mag- 
netic impurities. Excluding that feature, the calculated 
temperature dependence of the spin susceptibility resem- 
bles the observation. Fisk et al. []To| found an empir- 
ical correlation between the lattice distortion and the 
spin susceptibility, which is intriguing since \T measures 
the effective moment squared of an isolated magnetic 
ion. Our theory of this functional correlation provides 
an alternative explanation in which both properties are 
driven by the same condensation phenomenon. More- 
over, the theory predicts a low-temperature indirect gap 
of [a(0)Vb] 2 /M^ ~ 44 K, in good agreement with the ex- 
perimental activation energy of 42 K determined from 
transport measurement JTl| . 
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FIG. 2. Temperature dependence of static Pauli spin sus- 
ceptibility and its correlation with lattice deformation (inset) 
for gd = 2.38, gf = 1.55 for the same model as in Fig. 1. 
Squares are experimental results from Ref. 10 for CeaBi4Pt3. 

Since the exciton condensation is intimately connected 
with the structural transformation, one would expect cor- 
responding changes in elastic constants and phonon prop- 
erties. With the addition of the deformation terms in 
Eqs. (0) and (||), to first order in the external stress in 
Eq. (|7|), the dynamics of the pseudospin and the boson 
condensate is governed by a set of five coupled linear 
equations. The linear response deduced from them in- 
volves the collective mode in general [Q. The elastic 
constant is C + C con d, with the second being from the 
response of the condensate. At the high-frequency sound 
wave limit, where w»T (the linewidth of the collective 
mode, taken to be of the order of the indirect gap) but 
u) <C 2aV (the direct gap) , the response involves the 
collective mode and the dynamic elastic constant is 



C* co „d= ^a 2 N s E 2 Q x 

4£> 2 - 4A£>! + X 2 D + 2AXV Q 2 (D D 2 - D 2 ) 
(1 + W 2 D X ) 2 - V 2 D {mV 2 D 2 - 2A) 

where A = for short, and 



D, 



k 



(0) 



(k)]' 



/(£ k _) - /(£ k4 



2E- 



k(^ 



(10) 



(ii) 



The relative velocity change Av s /v s — 'St(C con d)/2C is 
plotted in the inset of Fig. 3 for our simple model. For 
the parameters of Fig. 1, the unit of relative change is 
N s El/CW = 0.26%. It is evident that its tempera- 
ture dependence reflects the second-order phase transi- 
tion and differs greatly from the typical anharmonic ef- 
fect. In the static or low-frequency limit (ui <C T), the 
elastic response is decoupled from the collective mode 
and is given by the quasi-particle contributions. There 
is a characteristic step-like drop M in the inverse static 
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elastic constant at T c , with the magnitude of the drop 
proportional to S§. The static elastic constant is quali- 
tatively different from the high-frequency one. 
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FIG. 3. Relative phonon energy change and intrinsic 
linewidth in units of the dimensionless factor N a S^/GW for 
a phonon of unrenormalized frequency uj/W = 0.2 + il0~ 4 
for the simple model described in the text. Inset: Relative 
velocity change for the case of lattice softening for ultrasound 
of u/W = 1(T 6 . 

In the study of high-frequency phonons, we have taken, 
for simplicity, the electron-phonon interaction to be given 
by the same deformation potential as above. A phonon, 
with energy greater than the direct gap (u > 2a Vq) but 
limited by the Debye frequency (<~ 0.1 eV), can break 
particle-hole pairs which contribute to the imaginary part 
of Ccond- This is the well-known Landau damping. The 
relative phonon linewidth contains a term — Q[C l con d]/2C'. 
Correspondingly, the relative phonon energy change is 
given by 9fi[C con d]/2C. Figure 3 presents the numerical 
results for such a high-energy phonon within our simple 
model. The resonance behavior in both the frequency 
change and damping occurs when the phonon frequency 
equals the temperature-dependent direct gap. Neutron 
scattering by phonons jl3| should be able to detect such 
resonances. Similar behavior in superconductors has 
been predicted |Q and observed 

We have also considered a two-dimensional square 
lattice with unit cell size d and a conduction band 
e(k) = — (cos k x d + cos k v d). For hybridization V(k) — 
iVf)Sh\(k x d/2) between a d-electron at the corner and an 
/ in between, the ground state is a fully gapped insu- 
lator, similar to our simple model presented above. For 
hybridization L(k) = iVosm(k x d) between d-f electrons 
at neighboring corners the exciton condensate contains 
nodes in the excitation spectrum. Nevertheless the nu- 
merical results for a(T) and consequently the lattice de- 
formation resemble those in Fig.l. Its interesting low- 



temperature transport properties are beyond the scope 
of this Letter. 

In summary, the additional consideration of the ro- 
tational symmetry of the hybridization between d and / 
electrons to the finite-temperature, mean-field solution of 
the periodic Anderson lattice using the slave-boson tech- 
nique reveals the correlated insulator state to be a p-wave 
exciton condensate as well as the existence of a concomi- 
tant structural transformation. Theoretical deduction of 
the temperature dependence of the spontaneous lattice 
distortion and of the Pauli spin susceptibility have some 
experimental support. In particular, the empirical corre- 
lation JhJ between the lattice distortion and the spin sus- 
ceptibility in the correlated insulator Ce3Bi4Pt3 can be 
explained in terms of the temperature dependence of ex- 
citon condensation. Moreover, intimate relation between 
the symmetry of the exciton condensate and the lattice 
distortion leads to theoretical predictions of (1) char- 
acteristic static elastic constant jump across the phase 
transition, (2) temperature dependence of the high fre- 
quency ultrasound velocity, and (3) resonance behavior 
of phonon frequency and damping as the gap parameter 
varies with temperature. All these consequences of the 
theory could be easily subjected to further experimental 
tests in mixed-valent correlated insulators. 
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